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i— i Two-point correlators represented by either perturbative or non-per- 

f^ turbative integral equations in Euclidean space are considered. In general, 

£™^ it is difficult to determine the analytic structure of arbitrary correlators 

Q, analytically. When relying on numerical methods to evaluate the analytic 

D structure, exact predictions of, e.g., branch point locations {i.e., the multi- 

4^ particle threshold) provide a useful check. These branch point locations can 

be derived by Cutkosky's cut rules. However, originally they were derived 

^^ in Minkowski space for propagators with real masses and are thus not 

^ directly applicable in Euclidean space and for propagators of a more general 

C\| form. Following similar considerations that led Karplus et al., Landau and 

Cutkosky more than 50 years ago to the mass summation formula that 
became known as Cutkosky's cut rules, we show how the position of branch 
points can be derived analytically in Euclidean space from propagators of 
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very general form. 
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1. Introduction and motivation 



H More than 50 years ago, several studies on the analytic structure of 

Feynman amplitudes were published [2 EJ [3] that shed light on the occur- 
rence of the spectral threshold. Here we apply this long-known technique 
to two-point functions at one-loop order in Euclidean space. Thereby very 
general analytic structures of the loop integrands are permitted which allows 
to employ a large class of propagators known in closed form for perturbative 
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calculations, as well as to extend this to non-perturbative equations which 
need to be solved self-consistently. In principle this is straightforward, but 
as we applied this technique to various problems in recent studies [U El E] 
and are using it in an ongoing non-perturbative study of the quark prop- 
agator, we will provide details of the technique here. The applicability of 
Cutkosky's cut rules in Euclidean space is discussed also in [7J using the in- 
verse Stieltjes transformation. The arguments presented here, on the other 
hand, are based directly on the loop integral and the analytic structure of 
the integrand. 

To start with, consider a two-point function which is expressed in Eu- 
clidean momentum space at one-loop order. In the perturbative case, such 
a two-point function might be written as 

jHp) = mpm-p))^ = / *» ^=»> (i) 

J (2tt) 4 g(p,q, cos 9) 

where d — 4 denotes the (Euclidean) dimension, q is some loop momentum 
and 9 the angle enclosed by the momenta p and q. The non-perturbative case 
will be addressed in Section |3j For convenience, the integrand is expressed 
as a fraction of two functions / and g, which provides access to the poles of 
the integrand via the zeros of g(p, q, cos #), assuming that the zeros of g lead 
to integrable singularities. Note that in general such integrals are divergent 
and need to be renormalized, e.g., via BPHZ. In the following we assume 
that eq. ([!]) is a renormalized expression. 

In Euclidean space, the momentum-space integral might be written 
in d-dimensional hyper-spherical coordinates (here we consider d = 4, a gen- 
eralization to arbitrary d is straightforward). Regardless of the dimension 
d, the number of integrals in equation reduces from d to 2. We thus 
integrate over the square of the inner momentum y = q 2 and z = cos 9. The 
expression reduces to 

JZ{x) = [ dyy I dzy/l-z* fi f' y ' Z \ , (2) 

Jo 7-i g{x,y,z) 

with x = p 2 and £ 2 is the UV cutoff. Note that the factor of (27r) 4 , as well 
as the constant contributions of the trivial integrals have been absorbed 
into the function /. When x E C instead of x E Kq", a thorough analysis 
of the analytic structure in the ^-complex plane allows for the prediction 
of non-analyticities of the result. The latter we will identify with branch 
points. 

When x is not on the positive real axis, the integration in y is no longer 
straightforward and it becomes necessary to deform the integration con- 
tour. As long as it is possible to integrate from to £ 2 without hitting any 
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Fig. 1. Branch cut for complex values of y as obtained from equation |3|. 



obstructions like branch cuts, the result is analytic. Non-analyticities in x 
arise, when this is prohibited by obstructions in the y-plane. These can be 
poles but also branch cuts induced by the angle integration. 

The latter arise when integrating over the (integrable) singularities of 
the angle integral. How this comes about can be illustrated in a very simple 
example. Consider the following integral for complex values of y, where 
a,beR, 
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With the particular (but arbitrary) choice of a = —1,6 = +1 we find that 
a branch cut arises for y G [—1, 1], see fig. [Tj This branch cut is present 
because z picks up the singular values of the integrand. This situation is of 
course much simpler than for physical examples, but it already features a 
branch cut which is induced by the z-integral and shows up in the complex 
?/-plane. Furthermore, we learn that we do not have to perform the angular 
integral in order to predict the branch cut structure in the y-plane. It is 
sufficient to find all complex values of y for which the integrand becomes 
singular. Thus, in view of equation S, the parametrization of the branch 
cut in the complex y-plane for given x G C reads 



c{x) = yo(x,z) 



=+i 



(4) 



where yo(x, z) are the zeros of the denominator for given complex x, para- 
metrized by z, i.e., yo is a solution w.r.t. y of 



g(x, y, z) = 0, x E C, x const., z G [—1, 1]. 



(5) 



Solving equation ^ yields all branch cuts in a parametrized form as given 
in Q, as well as all purely ^-dependent non-analyticities. In general there 
are both, poles and cuts. Non-analyticities in x arise when the restrictions 
imposed on the contour of y do not allow a proper deformation of the 
contour. These are the points we are after. To identify them we need to 



find the values of x, where such a contour deformation is not possible. This 
happens when an endpoint of a branch cut coincides with a pole and/or an 
endpoint of the solely y-dependent non-analyticities of the integrand. From 
this one can also see how the discontinuity at the branch cut develops: 
Depending on which side of the branch point a point x lies, the integration 
contour has necessarily a different shape and the value of the integral is not 
continuous across the branch cut; see fig. [2] for an example. 

We can summarize the procedure to identify the branch points as follows: 

• Step 1: Find all non-analyticities of the y-integrand. =>* Cuts and 
poles in the y-plane. 

• Step 2: Find the points x c where the endpoints of the induced cuts 
coincide with other non-analyticities in the ?/-plane which depend only 
on y ^ Candidates for branch points. 

• Step 3: Check if the ^-contour can be deformed properly for any x c 
from Step 2. If not, x c is a branch point x^. 

2. Correlator with real masses 

To confirm the results from the mass-summation formula of [U EJ [3] , let 
us consider a generic correlator with two masses mi E R + and rri2 G M + , 

■^(P) = (0(PM-P)U=4 = \S\ Ti i ( ? ,P L g) 2. ST ' (6) 

J (27t) 4 ((p - q) z + m{){q z + m^) 

where f(q,p — q) is a function constructed out of the Lorentz invariants. 
Switching to hyper-spherical coordinates we have 
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Jo J-i (x + y- l^fx^fyz + m\) (y + m|) ' 



with the usual identifications x = p 2 , y = q 2 and z — cos 9. All con- 
stants have been absorbed into f(x,y,z). According to the procedure out- 
lined above, we need to determine the zeros of the denominator to get 
parametrizations for the cuts (Step 1). In this case the solution can be 
easily worked out analytically and yields 

yi = -ml, (8) 

2/2(^5 z) — —m 2 — x + 2xz 2 — 2\i —m 2 xz 2 — x 2 z 2 + x 2 z 4 , 
2/3 (x, z) — —m\ — x + 2xz 2 + 2J—m\xz 2 — x 2 z 2 + x 2 z 4l . 
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Fig. 2. The analytic structure in the ^/-plane for mi = 0.5 and vri2 = 1. Slightly 
above/below the predicted branchpoint of xo = —(mi +m2) 2 = —2.25 the contour 
can still be drawn, while for x = xq there is no continuous deformation possible. 

According to Step 2 we solve y<i{x,z) — y\ and ys(x,z) = y\ for z = ±1 
with respect to x. This yields the following branch point candidates: 



— (mi ± m2) 2 . 



(9) 



The solution with the minus sign does not interfere with the integration 
contour (see fig. [2]) so that we are left with 



x b = -(mi + m 2 ) 2 . 



(10) 



This agrees with the result of [H [2 [3] if we use it for Euclidean momenta, 
i.e., if we plug in negative squares of the masses there: 



x ~ —(mi + m^) 2 = ( iy m\ + iym^ 



For an example with complex masses see [5] . 
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3. Fitted correlators and non-perturbative cases: An outlook 

The procedure outlined in the section above proved to be useful in per- 
turbative studies where the analytic continuation was performed numeri- 
cally. The branch point location (s) obtained from this analysis served as a 
check for the numerics. In [4], the case of complex conjugate masses was 
studied, while in [5l [6] complex conjugate masses as well as a propagator 
featuring a cut with a singular endpoint were considered. For the com- 
plex conjugate poles, all three resulting branch points have been found by 
this procedure. Note that one of them is time-like and thus corresponds to 
the physical multi-particle threshold, while two are complex. Due to the 
existence of the latter no Kallen-Lehmann representation is possible. 

The next step is a self-consistent determination of propagators and their 
analytic structure in the non-perturbative regime of QCD, e.g., as input for 




Fig. 3. The quark propagator Dyson- Schwinger equation. 

bound-state equations. Several studies of the non-perturbative properties of 
Green functions using Dyson- Schwinger equations (DSEs) have been pub- 
lished, see [SI [9] for recent examples. In an ongoing study we are looking 
into the analytic properties of the quark propagator DSE in Landau gauge. 
Since the structure of the r.h.s. of the equation depicted in fig. [3] is an in- 
tegral of a similar form as in eq. ([2]), the same procedure can be used. The 
additional complication of such a calculation is that the integral depends on 
the quantity on the l.h.s. and the equation has to be solved self-consistently. 
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